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SUMMARY 


A theoretical analysis is presented for the slideout dynamics of a vehicle 
equipped with a tricycle skid-type landing-gear system. These equations are 
reduced to three degrees of freedom, and the results of numerical calculations 
are compared to flight-test data obtained from the X-15 research vehicle. A 
comparison of these results shows that the three -degree -of -freedom equations can 
be used to adequately predict the slideout distance, the direction of lateral 
displacement, and the approximate lateral displacement. A numerical study shows 
that the velocity at which the aerodynamic influence on the vehicle becomes neg- 
ligible can be predicted by the three -degree -of -freedom equations. 


INTRODUCTION 


The successful recovery of a space vehicle is completed by a landing which 
assures the safety of the personnel and equipment onboard the vehicle. In addi- 
tion to absorbing the impact loads, the landing system of the space vehicle 
which performs a horizontal landing must provide stability and controllability 
during the subsequent runout. These characteristics, which are obtained on most 
conventional aircraft by a landing -gear system utilizing a three -point arrange- 
ment and equipped with brakes and steerable wheels, have been extensively treated 
by many investigators (refs. 1 to 4, for example). The inherent stability of the 
tricycle landing-gear system and the extensive experience gained from use on 
conventional aircraft make this type of system attractive for space vehicles. 

A tricycle system consisting of two main gear equipped with skids, placed 
symmetrically to the rear of the vehicle, and a nose gear, equipped with either 
a skid or a wheel, located on the forward section was studied in references 5 
to 9 for use on space vehicles. Landing impact has been investigated both ana- 
lytically (ref. 7 ) and experimentally (refs. 6, 8, and 9)* Slideout (runout) 
results obtained from flight and model tests (refs. 8 and 9) verify the direc- 
tional stability of the skid-type tricycle system during slideout; however, these 
results pertain to a particular vehicle and landing-gear configuration and cannot 
be applied directly to the design of other space vehicles. In addition, these 
tests did not evaluate controllability during the slideout. 



In order to provide a method for evaluating the slideout characteristics as 
well as the controllability of future vehicles, this paper presents a comprehen- 
sive theoretical investigation of the slideout dynamics of a vehicle equipped 
with a skid-type tricycle landing-gear system. Results of numerical calculations 
obtained on a digital computer are compared with flight-test results in order to 
assess the validity of the theoretical investigation. 


SYMBOLS AND COEFFICIENTS 


b 

C D 


wing span, ft 
drag coefficient. 
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qS 
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©h 


rate of change of drag coefficient with horizontal- tail deflection, 
dCp 

per radian 

° 5 h 


lift coefficient. 
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qS 


rate of change of lift coefficient with horizontal- tail deflection, 

55 ;, pe r radlan 
d6 h 

, , . , „ „ . . , Rolling moment 

rolling-moment coefficient, 

qSb 
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damping-in-roll derivative, — per radian 

°2V 


rate of change of rolling -moment coefficient with yawing angular- 

dC z 

velocity factor, — 7-, per radian 
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dC, 


effective dihedral derivative, , per radian 




rate of change of rolling -moment coefficient with rate of change of 

dc z 

angle -of- side slip factor, , per radian 
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rate of change of rolling-moment coefficient with respect to aileron 

d Ci 

deflection, - — , per radian 
o5 q 


rate of change of rolling-moment coefficient with respect to vertical- 


tail deflection. 


hCi 
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, per radian 


Cm 

Ctoq 

Si 


pitching-moment coefficient. 


Pitching moment 
qSS 


pitching-moment coefficient at zero angle of attack 
longitudinal- stability derivative, =^~, per radian 

rate of change of pitching-moment coefficient with rate of change of 

c) Cm 

angle -of -attack factor, —rz, per radian 

°2V 


rate of change of pitching -moment coefficient with angle of sideslip, 

dCr, 


r— per radian 

dp 


Cm, 


rate of change of pitching-moment coefficient with pitching angular- 


dq 


m 


velocity factor, — —> per radian 

2V 


Cm ; 


■5h 
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rate of change of pitching-moment coefficient with respect to 

dCm 

horizontal-tail deflection, , per radian 


yawing -moment coefficient, 
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Yawing moment 
qSb 


“'n 




rate of change of yawing-moment coefficient with rolling angular - 


dq 
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velocity factor, — — , per radian 

Sv 


dq 
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damping-in-yaw derivative, , per radian 

°2V 
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dc 
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directional- stability derivative , , per radian 


rate of change of yawing -moment coefficient with rate of change of 

dC n 

angle -of- side slip factor, — r— , per radian 

°2V 


rate of change of yawing-mcanent coefficient with respect to aileron 


deflection. 


dC, 

d6« 
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, per radian 


rate of change of yawing -moment coefficient with respect to vertical- 

da 


J n 


tail deflection, r — , per radian 

d5 v 


side-force coefficient. 


Side force 
qS 


rate of change of side-force coefficient with rolling angular - 

dCy 

velocity factor, , per radian 

rate of change of side-force coefficient with yawing angular -velocity 
dCy 

factor, -pg, per radian 




da 


side -force derivative, , per radian 

rate of change of side-force coefficient with rate of change of angl 

dCy 

of- side slip factor, —7—, per radian 

°2V 

rate of change of side-force coefficient with respect to aileron 

dCy 

deflection, , per radian 


rate of change of side -force coefficient with respect to vertical- 

dCy 

tail deflection, — 
dS Tr 


per radian 
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mean aerodynamic chord, ft 


d distance between center of gravity and main gear point of attach- 

ment measured in y-z plane parallel to z-axis, ft 

di,d 2 ,djj distance between center of gravity and left main skid, right main 

skid, and nose skid, respectively, measured in y-z plane parallel 
to z-axis, ft 

Ffo friction force on the landing gear, measured in the ground plane, lb 

I Fhi ,^p>*hw friction forces on the left main skid, right main skid, and nose 
skid, respectively, measured in the ground plane, lb 

F v landing -gear reaction normal to the ground plane, lb 

F V1 ,F V2 ,F V „ reactions of the left main gear, right main gear, and nose gear, 
respectively, normal to the ground plane, lb 

force measured parallel to x-, y-, and z-axes, respectively, lb 

hi, h 2 distance between point of attachment and left main skid and right 

main skid, respectively, measured in the y-z plane parallel to 
the z Q -axis, ft 

h]_,lv> distance between point of attachment and left main skid and right 

main skid, respectively, measured in the y 0 -z 0 plane parallel to 
the z G -axis, ft 

1 I x ,I y ,I z moments of inertia referred to the x-, y-, and z-axes, respectively, 
slug-ft2 

I xz product of inertia referred to the x- and z-axes, slug-ft^ 

distance between center of gravity and main gear point of attachment 
measured in x-z plane parallel to the x-axis, ft 

Zjjjj distance between center of gravity and nose gear point of attachment 

measured in x-z plane parallel to the x-axis, ft 

Z 2 , Zjj distance between center of gravity and left main skid, right main 

’ skid, and nose skid, respectively, measured in x-z plane parallel 

to the x-axis, ft 

moments about the x-, y-, and z-axes, respectively, ft-lb 
i m vehicle mass, slugs 

p,q,r rolling, pitching, and yawing angular velocities, respectively, 

measured about the x-, y-, and z-axes, radian/sec 
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v w 

v x> v yX 


V ; V 7 
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y G 
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dynamic pressure. 


|pV 2 , lb/sq ft 


wing area, sq ft 

distance between center of gravity and left main gear and right 
main gear point of attachment, respectively, measured in the 
x-z plane parallel to the y-axis, ft 

distance between center of gravity and left main skid and right 
main skid, respectively, measured in the x-z plane parallel to 
the y-axis, ft 

time, sec 

translational components of relative velocity vector parallel to 
x-, y-, and z-axes, respectively, ft/sec 

components of wind velocity parallel to x-, y-, and z-axes, 
respectively, ft/sec 

components of wind velocity relative to X- and Y-axes, 
respectively, ft/sec 

components of center-of -gravity translational velocity parallel to 
the x-, y-, and z-axes, respectively, ft/sec 

magnitude of vehicle relative velocity, ft/sec 

vehicle relative -velocity vector, ft/sec 

velocity vector of vehicle center of gravity relative to the ground, 
ft /sec 

magnitude of relative -velocity vector measured in the ground plane, 
ft/sec 

magnitude of wind velocity, ft/sec 

velocity vector of wind relative to the ground, ft/sec 

components of center-of -gravity velocity, relative to x-, y-, and 
z-axes, respectively, ft/sec 

components of center-of -gravity velocity parallel to the x 0 - 
and y Q -axes, respectively, ft/sec 

vehicle weight, lb 


X,Y,Z fixed-earth cartesian coordinate system 

x,y,z cartesian coordinates of body axis system 
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1 Subscripts: 


cartesian coordinates obtained by rotating X-, Y-, and Z-system 
through Euler angle 1 I/; distances measured parallel to these 
axes, ft 

angle of attack, radians 
angle of sideslip, radians 

angle measured between Vq and x D -axis, radians 
angle defined by equation (ll), radians 
aileron deflection, radians 

horizontal-tail (elevator) deflection, radians 

vertical-tail (rudder) deflection, radians 

change in vertical height of left and right main gear, 
respectively, ft 

length of time increment, sec 

angle between the x Q -axis and the friction force on the left main 
skid, right main skid, and nose skid, respectively, radians 

coefficient of friction for a skid 

coefficient of friction for the main skids and nose skid, 
respectively 

density of air, slugs/cu ft 

Euler yaw, pitch, and roll angle, respectively, radians 
angle between wind velocity vector and X-axis, radians 


i refers to body position or skid being considered 

n,n-l,n+l time interval being considered in numerical integration, n=0, 1, 

2 , 3 . . . 

aero aerodynamic inputs only 

gear landing gear only 

touchdown quantity measured at touchdown 

A dot above a quantity denotes the first derivative with respect to time. A 
double dot above a quantity denotes the second derivative with respect to time. 


ANALYSIS 


Slideout Equations for a Vehicle With 
Three Skid-Type Landing Gear 

The equations of motion are presented for the slideout of a vehicle with a 
tricycle landing -gear system equipped with skids. The vehicle geometry, tody- 
axis coordinate system, and directions of friction forces and gear reactions used 
in this analysis are shown in figure 1. The x-, y-, and z-tody axes, originating 







X 


at the center of gravity, are 
chosen so that the x-z plane is 
the vehicle plane of symmetry. 

The friction force F h lies in 
the ground plane and the gear 
reaction F v is normal to this 
plane. 

The transformation from a 
fixed-earth system of coordinate 
axes X, Y, and Z to the body 

axes x, y, and z is illus- 

trated in figure 2. A rotation 
about the Z-axis through the 
Euler yaw angle f results in 
the coordinate system x 0 , y 0 , 

z 0 . The x, y, z system is 

obtained by rotating the x 0 , 
y Q , z Q system, first through 
the Euler pitch angle 0, and 
then through the roll angle 9. 
Thus, the ground-plane components 
of the x-, y-, z-axes lie in the 
x 0 , y 0 , z 0 system. The trans- 



Z, z 0 


Figure 2.- Orientation of body axes with respect to 
earth axes. 


formation is dependent upon the order of rotation which, for this analysis, is 
6 , and qp. 


y 0 


Referring to figures 1 and 2, and resolving friction forces and gear reac- 
tions into components parallel to the body axes, the following equations, re- 
ferred to the body axes, are obtained for the contribution of the landing-gear 
forces to the total forces and moments on the vehicle: 


= - (f^ cos + F h2 cos Xg + F hlJ cos X^cos Q + ^F Vl 


+ F VN )sin 0 


^Vgear " <P sin X 1 + sin 9 sin <P cos x l) " F h 2 ( cos <P sin ^2 + sin 9 sin ^ cos * 2 ) 

- F^cos cp sin + sin 0 sin cp cos Xjjj - ( Fv l + F v 2 + Fv n) C0S Q sin ^ 

ZF Z g ear “ sin 9 cos 9 cos X 1 - sin 9 sin x l) ** F h 2 ( sin 9 cos V cos “ sin 9 sin ^ 2 ) 

- Fj 1 j(sin 0 cos cp cos - sin cp sin - ( F V]_ + F v 2 + F vn) cos 9 cos 9 


1 


9 


EM Xgear = Fh l[( S Ai cos 9 + d l sin 9) sin 0 cos - ^S Al sin cp - cos cpjsin \j_ 

- F h 2 [( SA 2 cos 'P " d 2 sin 9) sin 0 cos Xg - (s Ag sin cp + d 2 cos cpjsin Xgj 

+ F hu d N( cos 9 sin *-N + sin 0 sin 9 cos ^-n) + F Vl ( S A ^ cos cp + dj. sin cp^cos 9 

- F v 2 (Sa 2 cos cp - d 2 sin cp^cos 0 + Fy^d^- cos 0 sin cp 

£% g ear = ^lj^l sin 9 sin X 1 ~ ( l l sin 6 cos 9 + d l cos 0) cos X-j_ 

+ F h 2 ^2 sin 9 sin ^2 " { l 2 sirL Q cos 9 + d 2 cos 0 ) cos ^2 

“ F hJ*N sin 9 sin **N " ( Z N sin e COs 9 “ d N cos 0 ) cos ^ 

+ F v ^ dj sin 0 - Zj_ cos 0 cos cpj + F v ^ d 2 sin 0 - Z 2 cos 0 cos cp^ 


+ F VjJ ^d N sin 0 + Z N cos 0 cos 


9 ) 


^kgear ~ *1 COs ^ sin \l + (^1 sin 0 sin 9 “ S Ai cos ^cos X^_ 

j 

+ Fh 2 ^ 2 cos cp sin X 2 + (l^ sin 0 sin cp + cos 0^cos Xg 

-F^Zn^cos cp sin Xjj + sin 0 sin cp cos X N j 

+ sin 0 + cos 0 sin cpj - Fyg^S^ sin 0 - l 2 cos 0 sin cp^ 

- Fvjj^u cos 0 sin cp 


(Id) 


(le) 


(If) 


Aerodynamic forces and moments, referred to body axes, are of the form 
normally used in airplane stability analysis (ref. 10) and are written as: 


^xaero = ^ S ( C L Sln 0 “ C D C0S 0 ) 


ZF- 


Yaero 


= 


°Y 5a 5 a + C Y &v S v + C Yp P + ^(°Y r r + c Y p P + c YpP) 


(2a) 


(2b) 
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2F„ 


^aero 


= -qs(c L 


cos 9 + Op sin 9 


’) 


(2c) 


mx aero = ^[^S^a + C l§ y b v + C l p* 3 + w(% P + °*r r + C ^)] (2d) 

2%aero = + C^a + + ^(Cm q q + C^d)] (2e) 

2Mz aero = q^Cng^a + c n &v 5 v + c npP + 2v( Cn pP + C n r r + G npP) J ( 2f ) 

The total forces and moments on the vehicle during slideout are obtained by 
adding inertia and weight contributions to the sums of equations (l) and (2). 
Thus, the equilibrium equations referred to the vehicle body axes are: 


2 F x““ ln (^ + 1 w - rv) + qS^Cx, sin 9 - Cp cos ej + ( F Vl + F V g + FvjJsin 9 

-^Fh 1 cos X Y + Fh g cos \q + Ph N cos Xjj^cos 9 - W sin 6 = 0 ( 3 a ) 

SF y - -m(v + ru - pw) + i s [°Y 6a 5 a + C Ye ^S v + C Yp 3 + 7 ^(%r + C^p + 0^(3 

- (Fv ± + F V g + F V jj^cos 9 sin 9 - Fj^cos 9 sin X^_ + sin 9 sin cp cos X-^ 

- F h 2 ( COS sin X 2 + sin 9 sin *P cos ^e) " F hu( cos <P sin X N + sin 9 sin <P cos *n) 

+ W cos 9 sin cp = 0 ( 3^) 


£F Z = -m(w + pv - qu) - qs^ cos 9 + C D sin 9j - ( F v± + F V g + F v ^cos 0 cos cp 

‘ ( s: *- n 6 cos <P °os X-j^ - sin 9 sin Xj_j - Fj^sin 9 cos cp cos Xg - sin qp sin Xg^ 

- 9 cos 9 cos X N - sin 9 sin Xj^ + W cos 9 cos 9=0 ( 3c) 
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™x = -^P + + pq) + ( J y - + + C lp P 

+ ^(% p + C *r r + C *p^)] + F V i(s^ cos cp + dj sin cp^cos 6 

- f v 2 ( S cos 9-^2 sin c p) cos 0 + F v^% cos & sin 9 

+ Fh^J^S^ cos cp + d^_ sin cp^sin G cos s ^- n 9 " ^ cos <pjsin 

- F^g ( cos cp - d 2 sin cpjsin G cos \ 2 “ ( S A^ s:i ' n 9 + d 2 cos 9j sin 

+ F^d^cos cp sin Xjj + sin 0 sin cp cos * 0 

= -3y4 + (^z " ^x)p r + ^-xz( r2 “ P^) + qSdJ^Cm^&h + PmQ + 

+ CmpP + q + C^a)] + Pv 1 (d 1 sin 6 - l ± cos 0 cos cp) 

+ Fvg^d .2 sin 0 - Z 2 cos 0 cos cp) + Fy^djj sin 0 + Zjy cos 0 cos cp) 

+ Pjj^Tzj sin cp sin X^_ - sin 0 cos cp + d^ cos 0)cos X-jJ 
Z 2 sin cp sin X 2 - ^ 12 sin 0 cos cp + d 2 cos 0)cos XgiJ 

Z N sin cp sin X N - ^Z N sin 0 cos cp - d N cos 0)cos X N J = 0 


+ F v 


th 2 

- 


EMz = 


- X z* + I xz ( P - q r ) + (*x - 3 y)p ( l + Cn 5 a S a + C n &v 5 v + Cn^P 
+ ^(Cr^r + Cn p p + Cn^p)] + F Vl (s Al sin 0 + l ± cos 0 sin cp) 

- F^J^S, sin 0 - Z 2 cos 0 sin cp)- F Vn Z n cos 0 sin cp 

+ Fh 1 ^ 2 _ cos cp sin X^_ + ^Z^ sin 0 sin cp - cos djcos X-J 

+ F h 2 [^2 cos cp sin X 2 + (z 2 sin 0 sin cp + cos 0)cos X 2 J 

- F^Z n ^cos cp sin Xjj + sin 0 sin cp cos X N ) = 0 


(3d) 


(3e) 


(3f) 


12 



The skid, friction force is proportional to the gear reaction normal to the 
ground and is expressed as: 


F h = pF y 


( 4 ) 


For most applications the coefficient of friction p is considered to be 
constant ; however, it may vary with velocity and temperature as shown in refer- 
ences 7 and 11. Since the main skids on a reentry vehicle are expected to be 
constructed of the same material, the coefficients of friction for these skids 
are considered to be identical. The nose-gear skid may be made of a material 
with a coefficient of friction different from that for the main skids. Thus, the 
friction forces, in the form of equation ( 4 ), may be written as: 


= J 

~ ^ 1^2 > ( 5 ) 

F h w = I 


The dynamics of transformations between earth-fixed and body coordinates 
yield the following equations relating body and Euler rotational, rates (see 
ref. 12): 


p = cp - f sin 9 


q = 9 cos cp + i/ cos 9 sin cp 


( 6 ) 


r = ^ cos 9 cos cp - 9 sin cp 


The relative velocity and the components of relative velocity are related 
through the sideslip angle and the angle of attack. By definition. 


and 


P = sin -1 ^ 

—"1 w 
a = tan - 
u 


( 7 ) 
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As a result of equations (7), the relative-velocity components and their 
time derivatives are: 


v = V sin p 


u = V cos (3 cos a 


w = V cos p sin a 

. 

v = V sin p + Vp cos p 

u = V cos p cos a - V(p sin p cos a + a cos p sin a) 
w = V cos p sin a - V(p sin p sin a - a, cos p cos a) 


In general, the landing gear -will deflect and bend as the loads on the gear 
vary during landing and slideout. The rate of change of total vertical gear 
deflection can be related to the vertical velocity of the center of gravity as 
shown in appendix A. The results of appendix A are used to determine the gear 
deflection which is related to the gear loading. One method of determining the 
relationship between the total vertical deflection and gear loading is to conducl 
tests of the actual hardware. Thus, the gear loading which causes a specific 
deflection can be found. In addition, the moment arms from the center of gravity 
to the skids can be corrected as the gear deflect and bend. 

The vertical velocity z Q of the center of gravity is related to the 
vehicle geometry and rate of vertical deflection of the gear by equations (A3) 
and (a 4) and to the relative -velocity components by (see ref. 13) : 



z n = -u sin 0 + v sin cp cos 0 + w cos cp cos 0 


(9! 


The angle of attack is related to the vertical velocity of the center of 
gravity by 

a = 0 - 7 


do) 


where 7, referring to the vehicle attitude, is expressed as 


7 = 


sin 


-1 


z 


o 


V cos p 


(li: 
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Additional equations are required if the elasticity of the vehicle struc- 
ture, as well as that of the landing gear, is to he taken into account. These 
elastic conditions are not considered in this paper. 


Simplified Slideout Equations for a Vehicle With 
Three Skid-Type Landing Gear 

A practical approximation to the equations of motion presented in the pre- 
ceding section may he obtained from the assumptions that the vehicle has rigid 
gear during slideout and that all skids must remain in contact with the ground 
at all times. These assumptions are valid if (l) the vehicle is sliding on a 
flat surface, (2) the initial landing -impact oscillations have heen damped, and 
(3) gear deflections re stilting from all other inputs are negligible. 

The above assumptions imply that the roll and pitch rates must be zero; 
therefore, the roll and pitch attitudes remain constant. In addition, since it 
has been assumed that the vehicle is symmetrical about the x-z plane and that 
the geometry of the main gear is identical, the roll angle must be zero. These 
physical restraints on the vehicle attitude during slideout are stated as 

0 = 0 , 9=0 I 

and / (12) 

0 = constant, 9 = constant = 0 J 


The symmetry of the vehicle and the assumption of rigid gear also imply 

that 


and 



( 13 ) 


Since the center of gravity cannot have any vertical motion as a result of 
the previously stated assumptions, and if it is assumed that any existing wind 
"blows parallel to the ground plane, the velocity of the vehicle relative to the 
airstream during slideout will be parallel to the ground plane. Consequently, 
equation (10) becomes 

a = 0 = constant ) 

and > ( 14 ) 

a = 0 ) 
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The substitution of equations (5), (6), and (8) into the equilibrium 
equations (3) yields, upon the application of equations (12) to (l4) , the 
following equations for the slideout of a vehicle with a tricycle, rigid, 
type landing-gear system: 

ZF X = mjv cos (3 - (p + i })V sin pj cos 0 + ^pV 2 S^C D cos 9 - sin 0^ 

+ cos 9 cos X^ ” sin + ^p^ cos 9 cos X2 “ sin ^ 

5 ) ! 


+ cos 0 cos X^ - sin 0jF y ^ + ^ sin 0 = 0 


ZFy = mjv sin p + (p + i)V cos p] - |pV 2 S 

+ ^( c y/ cos 0 - c y/ sin 0 + °yJ) 


+ p 2 sin X 2 Fv 2 + Hjj sin X n F Vn = 0 


°*6 a 5 a + G Y5 v 5 v + °YpP 


+ p^_ sin XnF- 


l r v-, 


SF Z = m|v cos P - (p + tf)V sin pj sin 9 + cos 9 + Cp sin ©j 

+ ^p-^ sin 9 cos Xj_ + cos 0^Fy_^ + ^p^ sin 9 cos Xg + cos 0^F V ^ 
+ sin 9 cos Xjj + cos - W cos 0 = 0 


ZMx = ^Ix sin 0 + I xz cos 9jf + ipV^Sb 


Cz 5 a &a + Cz 6 v &v + C *p P 


^j(c lr i cos 0 - C Z J/ sin 0 + C 2 .p 


cos 0 


Jj i l^ s A 1 sin 0 cos ^1 + d l sin ^l) + S A- 

P-^Sa^ sin 9 cos ^2 ” ^1 sin ^2 ) + gos 0 


v 2 


+ Pj,d N sin X^F^ = 0 


skid- 


(15a) 


( 15 ^) 


(15c) 


(15<1) 


16 


IMy = |^I X - I z )sin 0 cos 0 + l xz |cos 2 0 - sin 2 ©^^ 

+ |pV 2 Sc + PniQ + + ) 

+ Jd-^sin 0 - (i-^ cos 0 cos X^j - Z-^cos 0 + ^ s ^ n ® cos ^l) 

+ ^d-^sin 0 - |i^ cos 0 cos K^j ~ ^( cos ® + ^1 s ^ n ® cos ^ 2 ) 

+ |djj ^sin 0 - |ijj cos 0 cos \ N j + Z^cos 0 + t-tjf sin 0 cos X^ 

SMz « -(l z cos 0 + I xz sin ejif + |pV 2 Sb[^Cn 5a 5 a + C n6 5 V + 

+ ^( c n r ^ cos e ~ Cnp^ sin 0 + C n 3 ) J 
+ sin X^ - S A ^ cos 0 cos X-^ + S A ^ sin 

+ sin X 2 + S^ cos 0 cos - S A sin ®j^Vg 




F v _ * 0 (l5e) 


N 


- sin x n f v n - 0 


d5f) 


The trigonometric relations for X^, X 2 , and Xjj are derived in 

appendix B. These relations, given by equations (Bl) , (B8) , and (B9) , yield 
[ the angles as functions of relative velocity, sideslip angle, and body yaw rate. 
Thus, the six simplified slideout equations ( 15 ) are in terms of six variables: 


V, (3, t' , Fy-, j F V g, and F- 


' V N* 


Simplified Slideout Equations for a 
Vehicle in a Low Pitch Attitude 

The equations of motions for slideout can he further simplified if the 
vehicle is in a low pitch attitude, that is, low angle of attack. Since the 
pitch angle is small, the following simplifications can be used: 

1. The trigonometric relations may be written as 


sin 0 « 0 
cos 0 « 1 


IT 



2. The z-component of relative velocity is small compared to the x- and 
y-components . Therefore, 

v « 0 
W as 0 

3 . Retaining only significant terms, equations (6) and their derivatives 
become for the slideout condition 

P as 0 p 0 

q = 0 q = 0 

• • • • 

r «* Tp r ijj 

If the above simplifications are applied to the basic slideout equations 
( 15 ), they become: 

ZF X = (m cos 3 j V - (mV sin 3^(3 - mrV sin 3 + Tp(Cj} - ^ 

+ (p x cos - 0^F Vl + (p x cos \ 2 - ejF V2 + (p N cos A n - e)F VN 

+ we = 0 ( 16 a) 


2F y = 

= (m 

l sin 3^V 

+ (mV 

cos f 

3 + mrV cos 3 

- f[v ♦ 

Cy 8 V + Cy p 

Y s v y p 



+ 

s(v - 

°YpP)] 

v 2 + 


sin 


+ (p]_ sin A 2 j 

F V2 



+ 

(m- n sln } 


= 0 






( 16 b) 

2F Z = 

ps 

2 

( C L + C D £ 

^V 2 + 

(hl0 

cos 

X 1 + 1 ) 


+ (p 1 0 cos \ 2 

+ ’■K 



+ 

(p^e cos 

a n + 1 

K 

- W 

= 0 




( 16 c) 
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ZM* = I xz r + ^[c l5a 6 a + C Z6 v 8 v + C^fB + ^(c z/ + C^V 2 
+ ^1 (^Ai 0 cos X-l + d x sin xj + SaJf^ 

- [^l( S A 2 0 cos A 2 - d l sin *2) + S A 2 F v 2 + (^N sin A n)% = 0 


ZMy = 


Ixzr 2 + ^(Ca^Sh + % + <V + VK 

+ ^1(0 ~ cos x x ) - 1-^1 + hi 0 cos ^i)J F Vl 

+ " hq cos ^2) ~ l l( d + ^l 0 cos ^2) F v 2 

+ |j%^9 - |ijj cos Xj^j + + M-n^ cos Xjjj 


F =0 


2M z = 


+ f^s/a + + V + ^(V + c^p)]? 2 

[“it 

[ni( l l sln X 2 + % oos *2 ) - VK ‘ (“h 


+ I Ut (Iq sin X^_ ~ ®A-j_ X^ + F ^i 


Z N Sln A n) F v n 


= 0 


( l6d) 


(l6e) 


(l6f) 


Equations (l 6), upon application of the trigonometric relations from 
equations (Bl), (B8) , and (BIO), can he solved, yielding equations for V, 0, 
and r. Because of the complexity of the resulting nonlinear differential 
equations, numerical methods are used to obtain a solution. The finite- 
difference technique described in reference lk was used. Initial conditions 
and equations for determining the slideout distance are given in appendix C. 


RESULTS AND DISCUSSION 


The slideout dynamics of a space vehicle equipped with a tricycle skid-type 
landing-gear system are presented in the basic equations ( 3) • These equations 
describe the general motion of the vehicle during slideout, including rolling, 
pitching, and yawing motions, variable aerodynamics, and variable coefficients 
of friction. Additional motion caused by structural flexibility of the landing 
gear is considered in appendix A. The basic equations are simplified to three 
degrees of freedom in equations (15) if the vehicle has rigid gear and if all 
skids remain in contact with the ground at all times during slideout . These 
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simplifying assumptions can be verified from model tests for vehicles of various 
configurations. The complexity of equations (15), as a result of nonlinearity 
and coupled motion, can be reduced if the vehicle is at a low angle of attack 
during slideout as presented in equations (16) . The general motion equations (3) 
can be used to predict the slideout characteristics as well as the stability and 
controllability of the vehicle during slideout. Equations (15) and (l6) have 
been derived by assuming that the vehicle is inherently stable, that is, it will 
not pitch or roll. These equations can be used to predict slideout character- 
istics and controllability. 

In order to ascertain the validity of the theoretical investigation, results 
of calculations performed for equations (16) are compared to experimental slide- 
out data obtained from the X-15 research vehicle. The X-15 is shown during 
slideout in figure 3. Numerical calculations were made on a digital computer, 
using the geometry and aerodynamic coefficients of the X-15 listed in the 
following table: 



Geometry 

b = 22.36 ft 

i N = 23.20 ft 

c = 10.30 ft 

S = 200 sq ft 

d x = 5.184 ft 

S Al = 5.10 ft 

d 2 = 5-184 ft 

= 5.10 ft 

<% = 3.583 ft 

p. 1 = 0.28 

l ± = 15-90 ft 

Pfo = °.°4 

l 2 = 15.90 ft 

9 = - 2 . 35 ° 

Aerodynamics 

Cp = 0.120 where = 0° 

= -1.320 per radian 

Cp~ = 0.105 per radian 
D h 

Cn r " ^np = - 1.350 per radian 

C L = -0.050 where 6^=0° 

Cnp = 0.280 per radian 

Cl_ = 0.975 per radian 
°h 

C n = 0.077 per radian 
°a 

0 Zr - = -0.020 per radian 

C n ^ = - 0.149 per radian 

C 2^ = 0.023 per radian 

^Y r “ C'Yp Considered negligible 

C? R = 0.097 per radian 
°a 

Cy^ = per radian 

= 0.032 per radian 

Cy~. = -0.100 per radian 
b a 

+ C^a = per radian 

C^ Considered negligible 

Cy^ = 0.240 per radian 
e v 


I 
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Figure 3.- Attitude of X— 1 5 airplane during slideout (photo taken from motion-picture sequence). 


Positive directions for the aero- 
dynamic coefficients are shown in 
figure k. Equations (l 6) are made 
more amenable to calculation pro- 
cesses by assuming that during 
slideout all aerodynamic derivatives 
and the coefficients of friction are 
constant. In addition, the products 
of inertia of the X-15 airplane are 
negligible compared to the moments of 
inertia. The trigonometric relations 
given in appendix B are valid for the 
main gear; however, a fully castering 
nosewheel which alines itself to the 
direction of relative motion of the 
center of gravity is used on the X-15. 
Thus, for the nose gear, the angle 
Xjj is approximately equal to the 
side slip angle (3 . 

Numerical calculations were made 
by utilizing slideout data from the 
X -15 research vehicle. A flight was 
selected in which the pilot varied the 
control positions during the slideout 
to determine the effect of control 
inputs on the stability and steering 
capabilities. The time histories of 
aileron, horizontal tail, and vertical- 
tail deflections from this flight are 



Figure 4.- Positive sense for coordinate system, 
aerodynamic coefficients, and control deflections. 
Arrows indicate positive sense. 
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shown in figure 5* Since the digital program could not accept continuous control 
changes, it was necessary to approximate the control inputs as shown hy the 
dashed lines in the figure. 





Figure 5.- Time histories of control deflections. 

Results of the calculations are compared with the flight data in figure 6 
as the variation of lateral displacement from the initial heading at touchdown 
with the slideout distance. The control inputs given in figure 5 produce the 
perturbations in the lateral displacement shown in figure 6. The numerals along 
the curves of figure 6 indicate the time during slideout by which the control 
inputs and lateral displacements can be correlated. Comparison of the calculated 
results and flight data indicates that the effects of control inputs on the 
lateral displacement during slideout can be predicted reasonably well by 
equations (l 6) . One heading change not predicted by the calculations occurred at 
an actual slideout time of approximately 20 seconds. However, in figure 5 it 
can be seen that there were no recorded control inputs to change the heading 
during this portion of the slideout. This displacement was probably caused by a 
variation in the lakebed surface on which the X-15 lands. 

In addition to the effects of control inputs on slideout, it can be seen 
from figure 6 that the calculated slideout distance is 6,863 feet, which is 
5 percent greater than the actual slideout distance, and the calculated lateral 
displacement is 117 feet, which is 1 9 percent less than the actual value. These 
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results show that the three -degree-of -freedom slideout analysis represented by 
equations (l6) adequately predicts the slideout distance, the direction of lat- 
eral deviations, and the magnitude of lateral displacement of the vehicle. 

Several factors omitted from consideration in the calculations were: 
variations in the aerodynamic derivatives and the coefficients of friction; the 
effect of ground plane on the aerodynamic derivatives; the oscillations of the 
vehicle during the landing impact; and deviations from the initial heading, 
although no intentional control inputs are made. A detailed study of the influ- 
ence of these factors on the slideout is beyond the scope of this paper. 

The effectiveness of lateral control inputs on the slideout is seen in the 
lateral displacement which is achieved. The results of a study of the effective- 
ness for one aileron setting for the X-15 are presented in figure 7* The lateral 
displacement is shown in the figure as a function of the velocity at which a 
control is removed during the slideout. The calculations were performed assuming 
that 15 0 of aileron control are applied immediately at a touchdown velocity of 
1 200 knots. The abscissa of figure 7 , therefore, is the velocity at which the 

control position is neutralized; that is, the point at zero indicates that the 
control is not removed during the slideout, whereas the point at 200 knots indi- 
cates that the control is removed immediately at touchdown. Figure 7 shows that 
if the control input is not removed until the vehicle has decelerated to about 
I 100 knots, there is negligible effect on the lateral displacement; however, if 
the control is removed immediately after touchdown, or before the vehicle has 
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decelerated to about 100 knots, the effect on the lateral displacement is signif- 
icant. These results indicate that the effect of control inputs is greatest 
between the touchdown velocity and 100 knots. Thus, the aerodynamic influence or 
the vehicle as indicated by the control effectiveness becomes negligible at about 
100 knots. 



Figure 7.- Lateral displacement as a function of the velocity at 
which an aileron deflection of 15° is removed. 

CONCLUDING REMARKS 


A theoretical analysis has been presented for the slideout dynamics of a 
vehicle equipped with a tricycle skid-type landing-gear system. The general 
slideout equations were reduced to three degrees of freedom, and calculations 
were made with these equations in order to compare the theoretical slideout with 
flight-test data and to study the effectiveness of aerodynamic controls on the 
lateral displacement during slideout. 

The results of the investigations indicate that the three -degree -of -freedom 
equations adequately predict slideout distance, the direction of lateral dis- 
placement, the approximate lateral displacement, and the velocity at which the 
aerodynamic influence on the vehicle becomes negligible. 


Flight Research Center, 

National Aeronautics and Space Administration, 
Edwards, Calif., March 1, 1963. 
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APPENDIX A 


DETERMINATION OF LANDING-GEAR DEFLECTION 


The vertical distances from the ground plane of the center of gravity and 
of the point of attachment of the main gear to the fuselage are related to the 
vehicle geometry as shown in figure 1. Analysis of the right main gear yields 


where 


= h r 


(d cos cp + S 2 sin cp^cos 0 + Ljj^ sin 0 


(Al) 


h 


2 


= hr 


cos 0 


Since the rate of change of height of the main gear point of attachment 
with respect to the ground is equal to the rate of gear vertical deflection, 
that is, 

h 2 = &m 2 (A2) 

the following relation is valid: 


= z Q - j^Ljj|yj cos d - (d cos cp + S 2 sin cp^sin 0 
- ^>2 cos cp - d sin cp^<j> cos 0 


(A3) 


The rate of vertical deflection for the left main gear is: 


= - jjkfflyi cos d - (d cos cp - S x sin cp^ sin 0 

+ ^Sj_ cos cp + d sin cpj< p cos 0 


(A4) 


25 



APPENDIX B 


TRIGONOMETRIC RELATIONS FOR THE ANGLE X 


The angle in the ground plane between the x G -axis and the friction force on 
a skid is denoted by X-^ where the subscript i indicates the skid being con- 
sidered (see fig. 1). The friction force on each skid will be in a direction 
parallel and opposite to the direction of the velocity of the skid relative to 
the ground. If the magnitude of the skid velocity is known, it is possible to 
determine the trigonometric relations for the Xj_ angles. 

The relative velocity of any position on a vehicle as measured in the 
ground plane is expressed as: 


K 

)\* 1 


V* 

x o 

and 


/ 

sin cp 

/ 

•H 

+ V 


(Bl) 


The ground velocity components V^ o and Vy„ are related to body axes by 

(< 0 )i = v ii cos 0 + ( v y± sln ' 

(Vy Q )i = Vy ± cos cp - V'. sin cp 


(B2) 


where 


v Xi = u ' + qzj_ - ry ± 
Vy ± = y' + rx t - pz t 
V'z ± = w' + py ± - qx i 


(B3) 


The components of the relative velocity of the center of gravity are denoted 
by the primed quantities. In order to determine these velocities, it is first 
necessary to investigate the relationship of the vehicle relative velocity, 
center-of -gravity relative velocity, and wind velocity, which is 

V = v' + v w (b4) 
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where 


V is the velocity vector of the vehicle center of gravity relative to the 
airstream 

V 7 is the velocity vector of the vehicle center of gravity relative to the 
ground 

Vw is the velocity vector of the wind relative to the ground 


Assume that a wind of constant velocity V w is blowing at a constant angle 
to the initial heading X and parallel to the ground plane. Then, the wind 
components in the fixed-earth system are: 


Uw q = V w cos * w 
v w Q = v w sin *w 



Transformation of 
u w = (%, cos * + 
v w = (v Vq cos i|r - 

w w = (% cos * + 


these velocities to body coordinates (see fig. 2) 
v Wq sin \(fjcos 0 

u^ sin \|r^cos cp + ^u^ cos \|^ + v w sin \|rjsin 0 sin 
v w ^ sin t^sin 9 cos cp - ^v w cos - u w sin \|rjsin 


yields: 
cp > (B 6) 

T I 


The components of wind velocity given by equations (B 6) and the components 

of V given by equations (8) are used in equation (B4) to determine the primed 
velocities as follows: 


u' = 


v = 


V cos p - V w ^cos i|f w cos \Jr + sin f w sin t^Jcos 9 

V sin p - V v [( 5in cos i|r - cos f w sin \(^cos cp 

+ (cos ijr w cos t + sin tw sin t)sin 9 sin <pj 


w' = V cos p sin 0 - V w 


(B7) 


(cos \|r w cos \|r + sin i|f w sin \|r^sin 0 cos cp 


^sin i)^ w cos i|r - cos sin 'I'jsin cp 
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The velocities given by equations (Bl) and (B2) are used to determine the 
trigonometric relations for X q, which are: 


cos X ^ = 


sin X^ = 


( v u 

( v °)i 

( v y°) 

( v «)i 


For the slideout of the vehicle represented by equations (15), the velocity 
relations given by equations (B2) become: 

= V cos 0 - t|ry ± - V w ^cos t w cos i|r + sin t w sin ) 

| (B9) 

( v yo) • = v sin ^ + ^( Xi cos 6 + z i sln 0 ) ” Vw ( sin cos ^ " cos sin *)) 

For the slideout of a vehicle in low pitch attitude, equations (B2) become: 

( V x 0 ) = V cos 0 - - V w ^cos \|r w cos \|r + sin t w sin I 

° 1 > (BIO] 

(v^q) = V sin 0 + ilrx i - V w ^sin t w cos \|r - cos t w sin \| r'j j 


28 



APPENDIX C 


INITIAL CONDITIONS AND CALCULATIONS OF EULER YAW ANGLE, 
SLIDEOUT DISTANCE, AND LATERAL DISTANCE 


The first derivative of V, (3, and r with respect to time can he written 
in the form of the parabolic -difference equations of reference 14- as 


^n+1 2e ^n + ^n-1 \ 

Pn+l = 2£ Pn + Pn-1 [ < Cl) 

r n+l “ 2e t> + r n-l ) 


where 


e is the length of the time increment 
n is the time interval being considered 


By designating the touchdown increment by n = 0, the conditions immediately 
after touchdown are dependent upon the conditions prior to and at touchdown. 

The conditions at touchdown are selected arbitrarily. It is assumed for the 
numerical example that the vehicle has a sideslip angle due only to a deviation 
in wind direction from the vehicle's initial heading and that the vehicle has no 
initial rotational rates. It is further assumed that these conditions hold 
immediately prior to touchdown. These conditions, together with the velocity 
initial condition, are expressed as: 

V 0 = ^touchdown \ 

e 0 = e_i^w [ ( C2 ) 

r 0 = r_i = ° / 


It is also assumed that the vehicle is decelerating linearly from a time 
prior to touchdown to a time after touchdown. Then V_-^ becomes: 


V -1 = 2V 0 - V 1 


(C3) 
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The equations for the Euler yaw angle, slideout distance, and lateral 
distance from initial heading, obtained by the use of the trapezoidal -difference 
techniques of reference l4 are: 


^n+1 = K + | ( r n+1 + r n) 


X n+1 = X n + | V Gn cos(^ n + p^) + V Gn+1 cos(^ n+1 + P^ +1 ) 


X n+1 = Y n + sin(* n + Pi) + V Gn+1 ain(* n+1 + P^ +1 ) 


The angle 0 is measured in the ground plane between the x Q -axis and the 
relative velocity of the center of gravity in the ground plane. The angle is 
determined by . , v 

= sln 'ty 

where referring to equations (12), (B2), and (B3), 


Vi = u f cos 0 + v r sin 0 


, ,/ t 

V v = v 
X °n 


and Vn is given by equation (Bl) . 
n 
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